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We investigate the thermopower of a metal-molecule-metal junction taking into account thermal 
effects on the junction. Based on analytical expressions and numerical simulations we show that 
the thermoelectric potential reveals valuable information on the mechanisms controlling the electron 
transfer process, including coherent transmission and thermalized hopping. We also show that at 
high temperatures the position of the Fermi energy relative to the molecular states can be easily 
deduced from the thermoelectric potential. Standard current-voltage measurements are insensitive 
to this information. 
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I. INTRODUCTION 

Understanding of charge transfer processes through 
sinde molecules is at the forefront of molecular electron- 
ics The electrical conductance of a single molecule 
coupled to metal electrodes has been recently measured 
by different techniques 0, . In these experiments the 
current-voltage (1-V) characteristic of the device is mea- 
sured, and information on the transport mechanism, the 
molecule-metal coupling strength, and the role played by 
nuclear motion in the conduction process are deduced. 

Experimental data ^ H, IS and theoretical studies 
suggest that two mechanisms are involved 
in electron transfer processes through molecular bridges: 
Super-exchange mechanism and thermal induced hop- 
ping. The super-exchange mechanism is a coherent (tun- 
neling) and short distance charge-transfer process, ft 
dominates transport when the molecular levels relevant 
for transport- Fermi level offset is greater than the ther- 
mal energy. In the opposite limit electron transmission 
occurs through sequential hopping along the bridge. This 
is a multistep process that allows long distance charge 
transfer. In a typical I-V experiment transport mecha- 
nisms are analyzed based on Arrhenius plot, revealing a 
characteristic transition from temperature (T) indepen- 
dent behavior at low T to a strong dependence at high 
T 0. 

It was recently suggested that thermoelectric mea- 
surements could provide additional insight into electron 
transport through single molecules, and contain informa- 
tion on the electronic and vibrational excitation spectrum 
of the molecule [11 [Hill. 

In a thermoelectric experiment the electric current in- 
duced by a finite temperature difference is investigated. 
Alternatively, the potential under the condition of zero 
current is measured. The thermoelectric power of mi- 
crostructures such as quantum dots , single electron 
transistors, [l^ and mesoscopic nanotubes [VA IT^ has 
been of interest since it yields directly the sign of the 
dominant charge carriers and the intrinsic conduction 
properties. In addition, it is sensitive to the electronic 
structure at the Fermi level. In two recent studies the 



thermoelectric voltage over a conjugated molecular con- 
ductor 12] and an atomic chain 13] was calculated. It 
yields valuable information- the location of the Fermi 
energy relative to the molecular levels. Inelastic inter- 
actions on the bridge were neglected in both cases. In 
a different work Koch et al. jl^] investigated the ther- 
mopower of a single molecule when both sequential tun- 
neling (lowest order tunneling process) and cotunneling 
(second order tunneling processes) take place. Vibra- 
tional features of the molecule were taken into account, 
but direct thermal activation of electrons onto the bridge 
and the possibility of diffusional transport were not in- 
cluded. In addition, this study was limited to the shortest 
molecular unit made of a single electronic state. 

In this paper we extend these ideas and analyze the 
thermoelectric voltage of a molecular junction of bridge 
length > 1 while taking into account both coherent 
and thermal interactions in a unified theory. Our model 
includes relaxation mechanisms on the bridge, arising 
from the interaction of the electronic system with other 
bridge or environmental degrees of freedom. The effect 
of these interactions is to open up a thermal channel for 
conduction. Hence in our model electrons can be trans- 
mitted through the bridge either coherently (tunneling), 
or by sequential hopping. We focus here on the non- 
resonant regime, where bridge energies are far above the 
chemical potentials of the metals. In this limit the actual 
population of electrons on the bridge is very small, and 
effects due to charging of the bridge (coulomb blockade) 
are negligible. 

We show through simple approximate expressions and 
with numerical examples that the functional behavior of 
the thermoelectric voltage is intrinsically different for the 
different modes of transfer, and that it clearly reflects the 
turnover between the two mechanisms. The thermoelec- 
tric phenomena can therefore serve as a significant tool 
for analyzing transport mechanisms complementing tra- 
ditional techniques. In addition, we find that when ther- 
mal interactions dominate the transport, the molecular 
levels -Fermi energy gap for transmission can be easily 
deduced from the thermopower value. This last property 
is crucial for interDretine: I-V results in transDort exoer- 
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iments. 



II. MODEL 

The model system consists of a metal-molecule-metal 
junction under an electrical bias and a thermal gradient. 
The molecule is described by a tight binding model with 
N sites with one state localized at each site. The first and 
last states are coupled to the left (L) and right (R) metal 
leads respectively. At equilibrium, ep specifies the Fermi 
energy of the two metals, taken as the zero of energy. 
Under a potential bias the leads are characterized by 
electrochemical potentials fj.^ and for the L and R 
sides. In addition, the two metals are kept at different 
temperatures Tl and Tr. For a schematic representation 
see Fig. 




FIG. 1: (Color online) Scheme of the model system: A 
molecule of A'^ sites connecting two metal leads, (a) No ap- 
plied bias, (b) e<j> > 0. 

We investigate the thermoelectric effect in this system 
using a variation of the model developed before to analyze 
thermal effects in electron transmission through molecu- 
lar bridges [^H^. There, the molecular system and the 
metals were in contact with a single thermal reservoir 
held at temperature T. Here, in contrast, we specify only 
the temperatures of the charge reservoirs (metals), and 
from this boundary condition the molecular temperature 
should be deduced. We explain below how to implement 
this modification into the formalism of Ref . . 

The system Hamiltonian includes five terms 

H = Hm + Hk + Hmk + Hb + Hmb- (1) 

Here Hm denotes the isolated molecule of length N with 
one electronic state at each site 

N N 

Hm = E +VY.{\J){J 1| + \j 1)01) • (2) 

i=i i=2 

The bridge energies Ej are taken to be equal at all sites at 
zero bias, E^ = E. V is, the nearest neighbors electronic 



coupling. The left and right charge carriers reservoirs 
(metals) are described by their set of electronic states 

HK = Y.^r\r){r\+Y,ei\l){l\. (3) 

reB. leL 

The leads are held at constant temperatures Tl (left) 
and Tr (right). In what follows we use the notations 
Ta = \Tl + Tr)/2 and AT = Tl - Tr. The molecule- 
metals interaction term is 

^A^^ =11 + I] K,jv|r)(^| -I- c.c. , (4) 

ieL reB. 

yielding the relaxation rate ri(e) = 2ttJ2i(=l ~ 
£;), and an analogous expression for Tr at the R side. 

The molecule is in contact with thermal degrees of free- 
dom, both internal (vibrations), and those related to the 
motion relative to the leads, all denoted as a "thermal 
bath" and included in Hb- The last term in the Hamil- 
tonian ^ describes the coupling of this thermal bath to 
the molecule. Specifically we use the following model for 
this interaction 

N 

ifMB-^i^,|j)01, (5) 

J = l 

where Fj are bath operators and \j) are molecular states, 
j — 1...N. The form of Eq. Q implies that the thermal 
bath induces energy dephasings on the local bridge sites. 
The bath operators are characterized by their time cor- 
relation function 

/oc />oo 
e^-HF,{t)Fy{0)) = / e-\F,,{Q)F,{t)), (6) 

where (3 — l/fcsT, T is the local thermal bath tempera- 
ture and ks is the Boltzmann constant. Assuming that 
thermal interactions on different molecular sites are not 
correlated, and going into the markovian limit, the cor- 
relation function becomes 

{F,{t)F,,{t'))^^5,,,,5{t^t'). (7) 

Here 7 is the dephasing rate reflecting the strength of the 
system-bath interaction: When 7 = 0, the electronic sys- 
tem and the thermal degrees of freedom are decoupled, 
and electrons move coherently along the wire. In con- 
trast, strong dephasing rates imply the dominance of the 
incoherent transmission mode [9j. It should be empha- 
sized that the bridge dephasing rate 7 and the electronic 
couplings between the sites V are in general temperature 
dependent. Since in the nonresonant regime thermal ac- 
tivation is the dominant temperature dependent factor, 
we simplify the discussion and ignore these corrections. 
The thermal bath temperature therefore enters the for- 
malism only through the detailed balance condition, Eq. 

The transport behavior of the system depends, among 
other factors, on the way the potential bias falls along the 
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molecular bridge. In what follows we use the following 
model for the electrostatic potential profile 

= CF + 60/2; M-R = - e(/)/2; 
El = E + 60/4; En = E- e(j)/4.; 
E, - E, j = 2..iV-l. (8) 

The qualitative aspects of the results presented below are 
not affected by this particular choice, since the relation 
Ai? = {E ~ €f) > e(j) is retained in this work. 

In the weak molecule-thermal bath coupling limit a 
computational scheme for evaluating the energy depen- 
dent transmission coefficient through a metal-molecule- 
metal junction for the Tl = Tr case was developed in 
Ref. Q. The method is based on the generalized quan- 
tum master equations extended to steady state situa- 
tions. Here we further extend this framework and study 
thermoelectric effects in molecular junctions. 



III. THERMOPOWER 

The current through the junction is calculated by gen- 
eralizing the Landauer formula 19] to situations involv- 
ing inelastic interactions on the bridge [20l l2l| 

Il^r = j d-'^o j d€fTL-,R{€o,ef)fL{eo) {1 - fnief)) , 

Ir^l = deo y"rfe/Tfl^i,(eo,e/)/i^(eo) (1 - /L(e/)) . 

(9) 

Here 7i^ij(eo,e/) denotes the transmission probability 
for an electron incoming from the left lead at the en- 
ergy eo to be emitted at the opposite side at ef. The 
energy difference ef — eo is exchanged with the thermal 
environment. This coefhcient depends on the molecu- 
lar parameters, the applied potential, the metal-molecule 
interaction, temperature and dephasing. It is not nec- 
essarily the same for the different directions. The cur- 
rent (defined positive when flowing left to right) is eval- 
uated by weighting the transmission probability by the 
appropriate combination of the Fermi distribution func- 
tions at the metals /i<-(e) = (e^^^""'"^) + where 
Pk = {kBTK)~^ is the inverse temperature {K — L,R). 
Unless specified, the integrals are all taken as J^^- 

The generalized Landauer equation does not take into 
account the effect of the contact population on the inelas- 
tic processes. Including such effects can be implemented 
by replacing the Fermi occupation factors by nonequi- 
librium electron distribution functions 0|- Yet it can 
be shown that Eq. @ well approximates the current 
when the transmission through the junction is signifi- 
cantly small (<< 1) m. Since we focus here on the out 



of resonance, small bias, weak electron-phonon interac- 
tion and weak metal-molecule coupling situation, trans- 
mission probabilities are always small and Eq. Q can 
be utilized. The inelastic Landauer formula can be also 
derived using the systematic nonequilibrium Green func- 
tion approach. Recent calculations of inelastic tunnel- 
ing currents yield an expression similar to © assuming 
weak coupling to the leads, weak electron-phonon cou- 
pling, and small bias 23|j. For recent discussions on the 
issue see Refs. (IJk |2^ |2j| . 
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FIG. 2; (Color online) Electron current vs. wire length 
iV for different temperatures: r=300 K (dashed), r=200 K 
(dashed-dotted) and T^lOO K (full). The other parameters 
are A£'=250 meV, V=25 meV, e(j>^2.5 meV, Fl = rjj=15 
meV, 7=1 meV. The inset presents the respective L to R 
transmission probability Tl^r for N=3 and eo ~ ep- r=300 
K ,200 K, 100 K top to bottom. 

We calculate next the current flowing through a bi- 
ased iV level molecular junction. We flrst consider the 
simple situation of = = T, i.e. we do not apply 
an external temperature gradient. Since this system was 
analyzed extensively in previous works, Refs. 0, l25j . 
we present here only the main results. Fig. [21 demon- 
strates the distance dependence of the current at different 
temperatures. The following set of parameters is used: 
AE = 250 meV, V = 25 meV, e0 = 2.5 meV and 7 = 1 
meV (Typical dephasing times in condensed phases are 
of order of 1 ps). The inset depicts the corresponding 
energy resolved transmission probability for an electron 
incoming at eo = £f for the N = 3 case. 

The following observations can be made: (i) The trans- 
mitted flux consists in general two main components: 
elastic tunneling at ef — eo and activated flux in the en- 
ergy range of the bridge states ^ AE. These two compo- 
nents can be clearly distinguished when bridge energies 
are high enough, AE ^ ksT and the coherent inter- 
action V is much smaller than AE. (ii) The coherent 
component is most important at low temperatures, large 
energy gaps and short chains. The incoherent contribu- 
tion dominates in the opposite limits, (iii) The electric 
current manifests a clear transition from an exponential 
decay typical to tunneling for short chains, to a weak (al- 
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gebraic) distance dependence for longer chains, charac- 
terizing thermal activation into the bridge. The turnover 
is shifted into higher N values when the temperature is 
lowered, (iv) The thermal component of the transmission 
is exponentially enhanced when increasing the tempera- 
ture. In contrast, the effect of the dephasing is more sub- 
tle. When the dephasing rate is very small 7 < 0.01 meV 
(using the parameters of Fig. [SJ, transport is coherent 
and the current increases with 7. For very high dephas- 
ing values, 7 ^ 100 meV, coherent effects are completely 
destroyed and the current goes down like I/7. See Ref. 
jPl for details. In this work we employ the intermediate 
values of 7 ^ 0.1 — 10 meV. 

Following these observations we can approximate the 
transmission coefficient by a generic functional form con- 
taining only its main features: For large energy gaps 
and for reasonable temperatures mixed coherent-inelastic 
contributions can be safely neglected and the transmis- 
sion is approximated by two separate terms 0, 



hopp ■ 



The first term stands for tunneling 



(10) 



(11) 



while the second contribution reproduces thermalized 
hopping through the junction 



Thoppieo,ef) = B(eo,e/)e 



-f3(EB-eo} 



(12) 



Eb is an effective bridge energy. It is related closely 
to the gap AE of the molecule in an equilibrium situ- 
ation, but it is modified by the molecule-metal interac- 
tions and the applied potential. The coefficients A and 
B depend on the molecule-metal coupling terms, the en- 
ergetics of the bridge, its length, and on the thermal pa- 
rameters T and 7. Here we assume that A does not 
depend on the thermal parameters, and that B is the 
same when exchanging between the initial and final en- 
ergies, B{eo,ef) — -B(e/,eo) ^ In what follows we uti- 
lize Eqs. (|11() - (|12|l for deriving simple expressions for the 
thermopower coefficient. 

We proceed to the relevant case of an electrically bi- 
ased junction under an additional temperature gradient. 
In our model the temperatures of the two leads L and R 
are kept fixed at Tl and Tr respectively, while the tem- 
perature of the molecular degrees of freedom is adjusted 
to the boundary conditions: At a steady state situation 
the temperature distribution along the molecule is deter- 
mined, among other factors, by the metals temperatures, 
the rate of energy dissipation on the molecule, and the 
rate at which energy is transferred away from the conduc- 
tor [2^ |2^ . Here we do not calculate this temperature 
gradient, but assume that at the left end of the molecular 
chain the temperature is nearly Tl, and similarly at the 
right edge it is close to T^. 

We assume next that the transmission Tk^k' depends 
on the temperature Tk (K=L,R), but is independent of 



Tfc'. Under this approximation the transmission coeffi- 
cient in Eq. ^ can be calculated using the procedure of 
Refs. [9ll25j without modifications, simply by employing 
the different temperatures when evaluating Tl^r and 
Tji^L. This approximation relies on the fact that in the 
K K' transmission process the temperature at the 
K end dominates the transport. This is true consider- 
ing both transport mechanisms, tunneling and sequen- 
tial hopping: Hopping through the bridge is triggered by 
thermal activation from the metal, see Eq. (|12|) . Assum- 
ing that this is the rate determining step, the relevant 
temperature is therefore Tk- The tunneling contribu- 
tion to transmission depends very weakly on the bridge's 
temperature 0. This naive approximation is adjusted 
by maintaining AT — > 0. Note that a full self consistent 
formalism (For example the Keldish- Kadanoff formalism 
[231 ) should naturally yield the temperature distribution 
on the junction without such assumptions. 

We define the thermoelectric voltage 0|/=o as the volt- 
age necessary for neutralizing the temperature induced 
current. The thermopower is the ratio of the potential 
energy e(j) to the temperature difference AT under the 
condition that the current vanishes 



S = — lim 



AT^O ksAT 



(13) 



We discuss next the limiting behavior of the ther- 
mopower ratio for the different transport mechanisms. 
The net current is calculated by considering separately 
the elastic (tunneling) and inelastic (hopping) compo- 
nents, Eqs. H10() - (|12|l . The tunneling current reduces to 
the standard Landauer formula ,19,] 



deo^(eo)[/L(eo)-/fl(eo)] 



(14) 



Assuming the transmission is a smooth function of the 
energy in comparison to ksTa, it can be expanded around 



^(e) = A{eF) + 
1 d^A 



dA 
Ik 



{e-ep) 



1 d^A 



2! 9e2 



3! de 



(15) 



In the linear response regime of small AT and cf) the Fermi 
functions in Eq. (|14|l are further expanded linearly 



/L(e)-/i?(e) 



df{e,eF,Ta) 
de 



e4>+ —(e-eF) 

^ a 



(16) 

Here the derivative of the Fermi function is calculated at 
the average values and cf- When we consider only the 
first two terms in Eq. (|15|) in conjunction with Eq. H16|l . 
we obtain the standard lowest order expression for elec- 
tron current due to both electric bias and temperature 
gradient [Illlllll 



Itunn = r^(e_F)(/i - 



'k%Ta dA 



de 



AT. (17) 
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The current through the device is zero when the potential 
difference is set to 



17=0 



TT^klTa d{ln{A)) 



3e 



de 



AT. 



(18) 



When higher order terms in A are necessary, we utilize 
the Sommerfeld expansion and get a power law series 
in Tq for the thermoelectric potential 



\i=o 



(19) 



We can further consider an explicit expression for A. If 
the bridge energies lies well above /i^ and fiR, a pertur- 
bative treatment leads to the superexchange result [s^l 



V 



2N 



{Eb - eo) 



2N 



(20) 



Here V is the coupling matrix element in the bridge and 
T=rK is the relaxation rate to the K metal. We sub- 
stitute this relation into Eq. (|19|l and get the tunneling 
contribution to the thermopower 



|/=o 



fcsAT 



E (2iV)" 



= 1,3,5.. 



(Eb - epY 



(21) 



This expression, though approximate, provides us with 
the important features of St'- the inverse dependence 
with energy gap, and its enhancement with size and tem- 
perature. 

When the gap {Eb — ep) is large and the bridge size 
N is long such that A is practically zero away from the 
bridge energies, i.e. A{e) cx d{e — Eb), the tunneling 
contribution from ep approaches zero. Electron trans- 
mission occurs then mainly via electrons in the tails of 
the Fermi distributions in the leads at energies around 
Eb- In this ballistic regime the current, Eq. H14|l . is zero 
when fL{EB) = fniEB) or 



(22) 



Using Tl = Ta + AT/2, Tr = Ta-AT/2, = ej. + e(/./2 
and /ifl = Ci? — 60/2, it reduces to 



' 7=0 



Ef 



kBAT 



knTa 



Sb- 



(23) 



This is the thermopower ratio in the ballistic regime: 
Electrons physically populate the molecule, but inelas- 
tic effects on the bridge are neglected. The thermopower 
in this case scales like T~^, in accordance with previous 
studies [T^. 

Next we estimate the thermopower when thermal in- 
teractions govern the transport across the bridge. We 
substitute Eq. into Eq. using B = Bk for the 



Tk->k' calculation, and get the net thermal (hopping) 
current 



'hopp 



def / deoBieo,ef)il - /L(eo))(l - /7?,(e/)) 



(24) 



It is zero when the term in the square parentheses van- 
ishes, i.e. 



Pl{Eb - A^l) = Pr{Eb - MTj), 



which leads to 



'7=0 



kBAT 



Eb — cf 
knTn 



= S 



H- 



(25) 



(26) 



This result is equivalent to Eq. H23|) . There we con- 
sidered resonant-coherent-transmission through a long 
chain, when tunneling from the Fermi energy is negligi- 
ble. Therefore, we cannot distinguish in the thermopower 
between transport due to thermal activation from the 
lead to the molecule and band motion of electrons at the 
tails of the Fermi function. We refer to both as thermal 
mechanisms. 

We compare next the tunneling thermopower term St, 
Eq. (|^ . to the thermahzed behavior Sh, (Eqs. if^ . 
(|26|l 'l. The following observations can be made: (i) St ~ 
{Ta/AEY, while Sh ^ AE/Ta- Here AE = Eg ~ ep- 
(ii) The tunneling term depends on the length of the 
molecule N. (iii) In both cases S does not depend on the 
molecule-metal interaction strength, given by the param- 
eter F. (iv) Measurement of S vs. T^ should yield the 
effective energy gap for transmission, and also hint on the 
transport mechanism. We expect that when increasing 
the temperature, the thermopower should first increase 
(tunneling behavior governs at low temperatures), then 
decay in a fashion when thermal activation domi- 
nates electron transfer. 



IV. 



RESULTS 



We investigate the thermopower behavior in our model 
system, Eqs. CP!-®! and show how the approximate ex- 
pressions, Eqs. (|21|) and H26|l . agree well with the nu- 
merical results. We also extract important energetic and 
dynamic information from the thermoelectric potential. 
The following set of parameters is used: N 1-5 units, 
molecular energies E ~ 100-500 meV (e^ is taken as 
zero). We focus on the limit E > V using V=25 meV. 
The metal-molecule relaxation rate is taken equal at the 
L and R sides, F^ = 15 meV. The system bath interac- 
tion is given by the dephasing parameter, 7 ~ 0.5 — 10 
meV. 

Fig. O shows the I-V characteristic for a chain of size 
N—A for different energy gaps E using Tl=335 K and 
rR=:300 K. Within this set of parameters the current is 
dominated by thermal effects, see Fig. Note that the 
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FIG. 3: I-V characteristics of the N=4 junction. The system 
parameters are 7=5 meV, Tl =335 K, rH=300 K. E = 315 
meV (dashed), E = 250 meV (full). The dotted line shows 
for reference the I—Q function. Inset: E — 185 meV. 




FIG. 4: The thermoelectric voltage plotted against AT/Ta 
for different metal-molecule energy gaps. J5=375 meV (full), 
£=315 meV (dashed), and £=250 meV (dashed-dotted). 
N^4, 7=5 meV, AT= 30 K. 
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FIG. 5: Dependence of the thermoelectric voltage on the 
bridge energy (E) in the superexcnahge regime. 7=0.5 meV, 
r„=100 K, Ar= -20 K. iV=2 (full), N=3 (dashed). Inset: 
iV=4. 



I-V characteristic is expected to be linear for e0 < ksTx 
for both coherent and incoherent modes of transfer |25j | . 
This property enables easy and accurate evaluation of the 
4>\i=o value by a linear fitting of a few I-V data points. 

We find that the results perfectly agree with Eq. H2t)|l . 
For _E=315 meV (dashed) the current vanishes at 
e(/)|/=o=-32 meV. When utihzing Eq. H26|l this number 
provides an effective energy gap of i?B=290 meV. This 
value is in agreement with the lowest diagonalized bridge 
energy of ~ i? - ^=315-25=290 meV. The same be- 
havior is obtained for iJ=250 meV (full) where we get 
e(/)|/=o=-25 meV, leading to Eb^221 meV. For ^=187 
meV (inset) the zero-current voltage is e0|/=o=-18 meV 
and the resulting gap is i?B=163 meV. We note that the 
effective energy gap Eb calculated (or measured) based 
on the thermoelectric effect is the real gap in the system, 
taking into account naturally and consistently the metal- 
molecule interaction, the thermal effects and the applied 
bias. 

Fig. ^demonstrates the temperature dependence of the 
thermoelectric voltage within the same range of param- 
eters. It provides another evidence that in the present 
case transport is dominated by thermal activation as 
4>i=o oc AT/Ta, in agreement with Eq. H2t)|l . The slopes 
for the E = 375,315,250 meV cases are 316, 267, 217 
meV respectively, producing the effective gaps. Note that 
for different thermoelectric voltages the bridge energies 
are slightly varied according to Eq. Q. For accurate 
results measurements should be done at AT — > 0. 

We turn now to the low temperature regime where tun- 
neling is expected to dominate charge transfer. Fig. |S1 
presents the potential 0|/=o for various chain lengths as 
a function of the energy gap E using Ta=100 K and 7=0.5 
meV. Indeed we find that the thermoelectric voltage de- 
creases with increasing gap in contrast to its behavior 
(in absolute values) in figures |3| and 01 In addition, the 
thermoelectric voltage is larger for longer bridges. These 
observations are consistent with Eq. H21|l . 

Next we show that the temperature induced turnover 
between coherent motion to incoherent transmission is 
reflected in the thermoelectric potential. Fig. El exempli- 
fies this behavior. We observe mainly three regimes in 
the main curve: (a) The potential is almost zero, (b) 
e(t>\i^o increases strongly with T^. (c) Above the thresh- 
old of Ta ~ 150 K the potential saturates, then goes 
down like cx AT/Ta- We can explain these results as 
follows: (a) At very low temperatures the thermoelec- 
tric potential is close to zero as tunneling transmission 
through the bridge is very small. Ballistic motion is not 
significant because of the low population of electrons at 
the metals around Eb- (b) When the temperature is in- 
creased electrons populate energies in the metals above 
the Fermi energy and tunneling becomes more proba- 
ble. Band motion and thermalized hopping also begin to 
contribute, (c) In this regime transport occurs through 
physical population of the bridge, and the thermopower 
behaves in accordance with Eq. H26|l . To conclude- while 
in region (a) transport is dominated by coherent interac- 
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tions, in region (c) it is induced by thermal effects. The 
intermediate area (b) presents a regime where coherent 
effects and thermal interactions mix. The inset displays 
the Arrhenius plot of Ln current vs. inverse temperature 
for a representative applied potential, = 6 meV (The 
results do not depend on the applied voltage for e0 up to 
~ 30 meV). The different regimes are marked according 
to the main plot. A clear transition at Ta ~ 100 K is ob- 
served: The current becomes temperature independent 
for lower temperatures since tunneling dominates. The 
behavior at region (c) is also clear: Ln(/) oc I/Tq. In con- 
trast, the behavior at the central area is obscure, and no 
clear transition is observed at Ta ~ 150 K. Thus, ther- 
mopower measurements may complement standard I-V 
studies, yielding valuable information about the junction 
energetics and charge transfer mechanisms. 
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FIG. 6: Temperature dependence of the thermoelectric volt- 
age. ^=250 meV, iV = 4, 7=1 meV, AT= -20 K. Inset: 
Arrhenius plot of current vs. inverse temperature for £</!> = 6 
meV. 
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FIG. 7: Temperature dependence of the thermoelectric volt- 
age for different bridge lengths. N=4: (full), A''=3 (dashed), 
iV=2 (dashed-dotted). £=250 meV, 7=1 meV, AT= -20 K. 



Finally, Fig. [7| displays the thermoelectric potential vs. 
temperature for molecules of different sizes: iV=4 (full), 
N=3 (dashed), and N=2 (dashed-dotted). We find that 
for short chains the transition points between the three 
regimes are shifted to higher temperatures as coherent ef- 
fects persist. For the N=2 chain transport is still largely 
controlled by tunneling at room temperature. Another 
important feature of the thermoelectric potential is its 
independence on TV at high temperatures, while it is 
strongly size dependent at low Ta- This observation is 
consistent with expressions (|^ and 



V. SUMMARY 

Using a simple model for electron transfer through a 
molecular junction, taking into account both coherent ef- 
fects and thermal interactions in the molecule, we have 
calculated the thcrmopower of the device assuming weak 
molecule-leads and weak molecule-thermal baths interac- 
tions. We have shown that the thermopower can provide 
information on charge transfer mechanisms, and that it 
can reveal the location of the Fermi energy relative to the 
molecular levels. Thermopower measurements can thus 
complement standard TV experiments, with the advan- 
tage of their insensitivity to the metal- molecule contact. 
For example, the thermopower can be used for compar- 
ing of transport experiments done on the same molecules 
but with different contact interactions and measurement 
methods f3L |. 
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